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■ Abstract 

'-^ . The study analyzes the temperature fluctuations in incompressible homoge- 

neous isotropic turbulence through the flnite scale Lyapunov analysis of the 
relative motion between two fluid particles. The analysis determines the tem- 
perature fluctuations through the Lyapunov theory of the local deformation, 
using the thermal energy equation. The study provides an explanation of the 
mechanism of temperature cascade, leads to the closure of the Corrsin equa- 
^ ! tion, and describes the statistics of the longitudinal temperature derivative. 

^ The results here obtained show that, in the case of self-similarity of velocity 

and temperature correlations, the temperature spectrum exhibits the scal- 
ing laws k", with n ^ —5/3, —1 and —17/3 -. — 11/3 depending upon the 
flow regime, and in agreement with the theoretical arguments of Obukhov- 
Corrsin and Batchelor and with the numerical simulations and experiments 
known from the literature. The longitudinal temperature derivative PDF is 
found to be a non-gaussian distribution function with null skewness, whose 
intermittency rises with the Taylor scale Peclet number. 

This study applies also to any passive scalar which exhibits diffusivity. 

Keywords: Lyapunov Analysis, Corrsin equation, von Karman-Howarth 
equation, Self-Similarity 



1. Introduction 

This work proposes the adoption of the flnite-scale Lyapunov theory, 
for studying the temperature fluctuations in incompressible homogeneous 
isotropic turbulence in an inflnite domain. The study is mainly motivated 
by the fact that, in isotropic turbulence, the temperature spectrum 9(k) 
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exhibits several scaling laws in the different wavelength ranges, depend- 
i ng on Tay l or scale Reynold s number R and Prandtl number Pr (^Corrsin 
fljAP 195lhJObukhovlfll949h . lBatchelOT] f ll959h . lB^tchelor et all f ll959h l This 
is due to the combined effect of these parameters which produces a peculiar 
connection between temperature fluctuations, fluid deformation and velocity 
field. 

For large values of R and Pr, ICorrsinI ( IJAP 195ll ) and lObukhovl (Il949( ) 
argumented, through the dimensional analysis, that Q(k) k~^I ^ in the so- 
called inertial-convective subrange (see Fig. [1]). iBatchelorl ( 1l959l ) considered 
the isotropic turbulence at high Prandtl number, when R is assigned. There, 
the author assumed that, at distances less than the Kolmogorov scale, the 
temperature fluctuations are mainly related to the strain rate associated to 
the smallest scales of the velocity field. As the result, he showed that ^ 
in the so-called viscous-convective interval, a region where the scales are less 
than the Kolmogoro v length (see Fig, [p . Dif fe rent experiments dealing wit h 
the grid turbulence f lCibson fc Schwar3in963l ) . iMvdlarski fc Warhaftl fll998l )) 
and calculations of the temperatue spectrum through numerical simulations 
( iDonzis et all (120 lOf ) and references therein) indicate that 6(k) follows the 
previous scaling laws. 

On the contrary, when Pr is very small, the high fluid conduc tivity deter- 
mines quite different situations with respect to the previous ones. iBatchelor et al 
(1l959l ) analyzed the small-scale variations of temperature fluctu ations in the 
case o f large conductivity, and found that Q{k) whereas iRogallo et al 

( 119891 ) determined the temperature spectrum through numerical simulations 
of a passive scalar convected by a velocity field with zero correlation time. 



Rogallo et all ( 119891 ) showed that, when the kinetic energy spectrum follows 
the Kolmogorov law E{k) ^ the temperature spectrum varies accord- 

ing to Q{k) k", with n ^ —11/3. 

Furthermore, experiments of grid turbulence state that both tempera- 
ture and velocity correlations are linked each other when Pr = 0(1) and 
that decay rate and characte ristic scales depend on the initial conditions. 
Specifically, iMills et all ( 1l958l ) obtained very important data about the air 
turbulence behind a heated grid. They carried out several measurements 
of nearly isotropic fluctuations of velocity and temperature at different dis- 
tances from the grid, and recognized that the temperature correlation fg is 
roughly equal to the longitudinal velocity correlation /, and that the triple 
correlation tempera ture- velocity is o f the o rder of the triple velocity cor- 
relation k. Later, IWarhaft fc Lumleyl ( 119781 ) experimentally showed that 
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spectrum shape and decay rate depend upon the initial conditions and that 
the mechanical-thermal time scale ratio tends to a value close to unity. 

Another important characteristics of is the self-similarity. This is 
related to the idea that the combined effect of temperature and kinetic energy 
cascade in conjunction with conductivity and viscosity, makes the tempera- 
ture correlation similar in the time. This pro perty was theoretically studied 



by George (see iGeorgd (Il988l ) , iGeorgd (119921 ) and references therein) which 



showed that the decaying isotropic turbulence reaches the self-similarity, 
where 6(k) is scaled by the Taylo r microscale whose c urrent value depends 
on the initial condition. Recently, lAntonia et all (120041 ) studied the temper- 
ature structure functions in decaying homogeneous isotropic turbulence and 
found that the standard deviation of the temperature, like the turbulent ki- 
netic energy, follows approximately the similarity over a wide range of length 
scales. There, the authors used this approximate similarity to calculate the 
third-order correlations and found satisfactory agreement between measured 
and calculated functions. 



From a theoretical point of view, the properties of G(k) can be investi- 
gated through the evolution equation of the temperature spectrum. Q{k) 
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i s the Fourier- Trans form of fg which in turn satisfies the Corrsin equation 
( ICorrsiru (jJAS 195ll )). This latter includes the term G, a quantity providing 
the temperature cascade, directly related to the triple correlation p*, and 
responsible for the thermal energy distribution at the several wavelengths. 
As G depends also on the velocity fluctuations, the Corrsin equation requires 
the knowledge of the velocity correlation /, thus it must be solved together 
to the von Karman-Howa rth equation. On this ar gument, some work has 
been written. For instance, Baev fc Chernykh ( 2010 i) (and references therein) 
studied temperature and kinetic energy spectra adopting a closure model for 
the Corrsin and von Karman-Howarth equations based on the gradient hy- 
pothesis, which incorporates empirical constants. Nevertheless, to the author 
knowledge, the estimation of Q{k,) based on the theoretical analysis of the 
closure of von Karman-Howarth and Corrsin equations has not received due 
attention. 

This is the motivation of the present work, whose main objective is to 
propose the closure of the Corrsin equation and a description of the statistics 
of the temperature derivative. The present study is ba sed on the finite-s cale 
Lyapunov theory, just used by Ide Divitiid (120 lOf ) and Ide DivitiisI (1201 if ) for 
determining the closure of the von Karman-Howarth equation and the statis- 
tics of the velocity difference. Here, this theory gives G in function of / and 
dfo/dr, and describes also the statistics of the temperature gradient through 
the Lyapunov analysis of the local strain and the canonical decomposition 
of temperature and velocity in terms of proper stochastic variables. The 
closure of the Corrsin equation is obtained considering that G is frame in- 
variant, thus G is calculated in the finite scale Lyapunov basis. The adoption 
of this basis is revealed to be an usefull choice for determining the analyti- 
cal expression of G. For what conc erns the vo i i Karm an-Howarth equation, 
the analytical closure proposed by Ide DivitiisI (120 lOf ) is here adopted. From 
the system of equations of von Karman-Howarth and Corrsin, an ordinary 
differential system is determined, through the hypothesis of self-similarity 
for / and fe. This differential system is first reduced to a Cauchy's initial 
condition problem, then it is numerically solved for several values of R and 
Pr. The results show that the temperature spectrum exhibits scaling laws 
whose exponents depends on R and Pr, in agreement wit h the exper i menta l 
and theoretic a l data of the litera t ure (ICorrsinI (I JAP 19511). Obukhovl ( 119491 ) . 
Rogallo et all (Il989h . iMills et all (Il958h . ICibson fc Schwarzl ( Il963h l As far 



as the statistics of the longitudinal temperature gradient is concerned, it is 
represented by non-gaussian PDF with null skewness and a Kurtosis greater 
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than three whose value rises with the Peclet number Pr R. 



2. Temperature correlation equation 

For sake of convenience, the procedure to obtain the Corrsin equation is 
here renewed. 

The isotropic homogeneous velocity and temperature fields are consid- 
ered, where fiuid viscosity u and thermal conductivity k are assigned quan- 
tities. The equations of the temperature fiuctuation ■(9 in two ponits x = 
(x, y, z) and x' = x + r, are 

dt dxk dxkdxk 



being r = {rx,ry,rz) the separation vector, x = k/{pCp) is the fiuid thermal 
diffusivity, and Cp is the specific heat at constant pressure. In case of ho- 
mogeneous temperature fiuctuations, the temperature correlation is defined 

as 

/. - (3) 



where 6 = is the standard deviation of the temperature fiuctuations, 

constant in the space. 

As well known, the evolution equation of fg is determined multipl ying Eg. 



([ID and (El) by {}' and {}, respectively, and summing the equations (I Corrsin 



( jjAS I951I )). The so obtained equation, averaged with respect to the ensam- 



ble of the temperature fiuctuations, leads to 

The first two terms of Eq. (|4]) express the time variations of fg and 6, 
the third one, arising from the convective terms, provides the mechanism of 
temperature cascade, whereas the last one, the laplacian in the spherical co- 
ordinates of fg, describes effects of the thermal diffusion. Because of isotropy 
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of temperature and velocity fluctuations, the third term can be expressed 
through the scalar G{r), an even function of |r| 



-G 



According to ICorrsinI ( 1.TAS 195ll ). ICorrsinI ( 1.TAP 195ll ). G is 



(5) 



(6) 



being p*(r) the triple correlation between temperature fluctuations in x and 
x', and the velocity component in x along the direction r 



p*{r) 



(7) 



and Ur is the longitudinal component of the velocity fluctuation. 

Now, the temperature cascade, responsible for the distribution of thermal 
energy at the different wave-lengths, does not modify 6', hence G(0) = and 
G fa near the origin. Therefore, when r — »■ 0, and for r = Eq. 

dlj) gives the evolution equation for 9 (ICorrsinI (jJAP 195ll )) 



dt 



(8) 



\g is the sc a le of the te mperature correlation, or Corrsin microscale, defined 
as f jCorrsinl fl.TAS 195 ih ) 



(9) 



where the superscript apex denotes the differentiation with respect to r. As 
the consequence, the evolution equation of fe is 



X 



G-2x 

whose boundary conditions are 

feiO) = 1, 

lim fg{r) = 



d'fe ^ 2dfe 

2 /y ^^'y 







(10) 



(11) 
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Note that Eq. ( ITU]) depends also on the velocity fluctuations through G(r) 
whose analytical expression is not given at this stage of this analysis. There- 
fore, Eq. ffTOj) is not closed and provides only a link between G and fg. 
Accordingly is related to /. As k and v are both constant, the von 
Karman-Howarth equation (see Eq. fl90|) . Appendix) is independent from 
fe and 6, hence Eq. flUUl) can be first solved separately, whereas Eq. f lTUl) 
requires the knowledge of / by means of G. 



3. Lyapunov Analysis of the temperature cascade 



The purpose of this section is to analyse the flow of the thermal e nergy 
cascade with the finite-scale Lyapunov theory used by Ide DivitiisI (120101 ) , and 
to propose an analytical expression for G which provides the closure of the 
Corrsin equation. To this end, consider now the expression of G (Eq. ([5])) 



G{r) 



d 



dri 



(12) 



This is frame invariant, therefore, for sake of convenience, G is expressed in 
the finite scale Lyapunov basis E\. This basis is associated to the problem 
of the relative motion between two fluid particles ( de Divitiis ( 2010l )) 



dp 

It 



u(x + p,t) - u(x, t), 
u(x,t) 



(13) 



where p gives the relative position between the particles, and u varies accord- 
ing to the Navier-Stokes equations. E\ is defined by means of the solutions 
pi, p2 and P3 of Eq. ( !T3|) . whose initial conditions Pi(0), P2(0) and ^3(0) are 
mutually orthogonal vectors which satisfy |pi(0)| = |P2(0)I = |p3(0)l = 
Specifically, Ex is obtained through the Gram-Schmidt orthonormalization 
process applied to Pi(t), p2if) and p^^it). 

In E\^ the velocity difference fluctuation is expressed as ( de Divitiis ( 2010l )) 



Au 



u 



u 



A(r)r + lja X r + C 



(14) 



where A is the maximal finite scale Lyapunov exponent (associated to the 
length r), defined by A(r) ^ 1/T dr/dt-r/r'^ dt, and calculated in function 
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Figure 2: Scheme of Finite-scale Lyapunov basis embedded into a material cylinder, at a 
given time 



of / as fide DivitiisI fl2010h ) 



A(r) = -v/2(W) 



(15) 



is the angular velocity of Ex with respect to the inertial frame of reference 
3?, and C = (Ci5C2,C3); related to the ot her two ex p onent s, makes Au a 
solenoidal field, and is expressed in E\ as ((de DivitiisI toiij} ) 

Cl = (Al-A)f?i, C2=(A2-A)f?2, C3=(A3-A)f?3 (16) 

where Ai, A2 and A3 are the Lyapunov exponents associated to the directions 
Qi, Q2 and ^3, respectively. As the consequence, A is a deterministic quantity, 
whereas uj\ is a fluctuating variable related to the relative motion between 
E\ and 3?. Without lack of generality, the coordinate qi is supposed to be 
associated to the maximal exponent, then Ai — )■ A, A2 = A3 = A^, qi diverges 
being >>> |^2|, 1^131, thus Ci Ci 7^ 



Ci = Ci, C2 = (Ac - A)^2, Cs = (Ac - A)^3 



(17) 



The exponents A2 = A3 = A^ are determined with the continuity equation. 
With reference to Fig. [21 the equation is written considering, at a given 
time, the mass balance associated to a material circular cylinder whose axis 
is parallel to the direction qi 



a Qi 
2- + — 
a Qi 







(18) 
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where the dot denotes the differentiation with respect to t, a is evaluated at 
the coordinate h Qi, with h G (0, 1), whereas Qi and a = \/ Q2 + oi are length 
and diameter of the cylinder. Therefore, Qi/ Qi and a /a identify A and A^, 
respectively, being 

Ac = -\ (19) 

Substituting Eq. into Eq. ( fT2l) . and taking into account that A is 

constant with respect to the operation of statistical average (i.e. (A..) = 
A(..)), G is the sum of three addends 



G = Gi + G2 + G3 (20) 



where 



Gi = -^m')^Qk), 
OQk 

G2 = -^m')Ck), (21) 

OQk 
d 

being Skij the Levi-Civita tensor. Gi reads as 

G. = -(,= (f Ar+^|(r'A)) (22) 
and G2 is written taking into account that A(r)r + ujx x r + ^ is solenoidal 

^2--^a + ^^(r A) (23) 

where r"^ = Qi + Q2 + Qs- For what concerns G3, it gives null contribution, as 
in isotropic turbulence, {dd'ui) is a function of r = |r| and this implies that 
G3 = 0. Therefore G is 

G . -A - f!?la (24) 
or OQk 
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The first term of Eq. fl24p is tlie consequence of A > and provides partially 
the mechanism of temperature cascade, giving a flow from small to big scales. 
The second one, related to the other two Lyapunov exponents, goes against 
the previous mechanism of cascade, being Au solenoidal, thus it will have 
the opposite sign of the first term. To obtain the second term, observe that 
into Eqs. f lT7|) . A2 = A3 = = —A/2, and this leads to 

C2 = o^Afc^ — , — Cs = o^Afc^ — [Zb) 

0()2 Z or r OQs 2 or r 

Moreover, because of the isotropy d{&^') /dgi(i must be of the form 

f^O ^ If A«'^^ (26) 
OQi 2 or r 

Hence 

Sm^^^_3,.dU,^ (27) 



drk 2 dr 

Accordingly, the analytical expression of G is 



G(r) = ^^Ar = e^uJ'-^ ^ (28) 
^ ' 2 dr y 2 dr ^ ^ 

Equation fl28|) . representing the thermal energy cascade, gives the proposed 
closure of the Corrsin equation and expresses the combined effect of tem- 
perature and velocity correlations. Its main asset with respect to the other 
models is that it is not based on the phenomenological assumption, but is 
derived from a specific finite-scale Lyapunov analysis, under the assumption 
of homogeneous isotropic turbulence. Equation (!28l) preserves 6 and provides 
the mechanism of the thermal energy transfer. This latter consists of a flow 
of the thermal energy from large to small scales which only redistributes the 
thermal energy between wavelengths and whose effe ctiveness d e pends upon 



/ and fe. According to the Lyapunov analysis of Ide Divitiid ( I2OIOI ). this 
mechanism can be viewed in the following manner. If, at an initial time to, 
a toroidal material volume E(io) is considered, which includes an assigned 
amount of thermal energy, its geometry and position change according to the 
fluid motion, and its dimensions will vary to preserve the volume. Choos- 
ing S in such a way that the maximal dimension of the toroid R increases 
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with t, the finite-scale Lyapunov analysis leads to i? ~ R{to)e'^^^~^°\ The 
thermal energy, initially enclosed into S(to), at the end of the fiuctuation is 
contained into whose dimensions are changed with respect to the initial 
time to- Therefore, the thermal energy is transferred from large to small 
scales, resulting enclosed in a more thin toroid. 



4. Formulation of the problem 

At this stage of the analysis, the problem for determining / and fe is 
formulated through the von Karman-Howarth and Corrsin equations, which 
are here reported 

Bf_m ^,J0^^4af\10.^ (29) 



dt \ dr"^ r dr J A| 

9fe_G{r) fd^fe 2dfe\ 12x 

whose boundary conditions are 

/(O) = 1, hm /(r) = 0, 



/,(0) = 1, hm fe{r) = 



(31) 



K and G are expressed by means of the finite scale Lyapunov analysis 
(Eqs. (1931) and (|28l)) 



3. /I-/ df 



or 



(32) 



ry q2 ~ f ^/f 



As k and u are considered to be assigned quantities, according to Eqs. 
fl29|) and fl30|) . fo is related to /, whereas / does not depend upon fg. 
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The energy spectrum E{k,) and th e transfer f u nctio n T{k) are the Fourier 
Transforms of / and K, respectively ( iBatchelon ( 1l953l )) 





1 r°° 








K{r) 



2 2 



sm Kr 



Kr 



cos Kr I dr 



(33) 



accordingly, the temperature spectr um Q{k ] and the temperature transfer 
function T{k) are here calculated as ( lOgural ( 119581 )) 



■ e{K) - 


2 












Kr sin Kr 


T{k) 




G{r) 





in such a way that 



/,(r) = / e(K)^^ dK, G(r) 



T{k)^ dK 

Kr 



and 



/•oo /"OO 

/ 0(/t) dK = ^^ / r(K) dK = o 

Jo Jo 



(34) 



(35) 



(36) 



5. Lyapunov analysis of the temperature fluctuations 

The proposed procedure for calculating the temperature fluc tuations is 
based on the Lyapunov analysis of the fluid strain just proposed by lde Divitiis 
( I2OIOI ). and on the adoption of Eq. ([1]). 

In order to obtain the temperature fluctuation, consider now the rela- 
tive motion between two contiguous particles, expressed by the infinitesimal 
separation vector d x which obeys to the equation 



die = Vu (ix 



(37) 



where (ix varies according to the velocity gradient wh i ch in turn follows the 
Navier-Stokes equations. As observed by Ide DivitiisI (I2OIOI ). in turbulence, 
(ix is much faster than the fluid state variables, and the Lyapunov analysis of 
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Eq. fl37p provides the expression of the local deformation in terms of maximal 
Lyapunov exponent A = A(0) > 

dx 



Now, the map x '■ xq — x, is the function which determines the current 
position X of a fl uid particle located at the referential position xq at t = to 
(jTruesdelll (119771 )). Equation ([T]) can be written in terms of the referential 
position xo 



dt \ dxop dxopdxoq dxh J dxh 

The adoption of the referential coordinates allows to factorize the temper- 
ature fluctuation and to express it in Lyapunov exponential form of the lo- 
cal deformation. As this deformation is assumed to be much more rapid 
than d'd/dxopUh and xd'^^/dxopdxoq, the temperature fluctuation can be ob- 
tained integrating Eq. fl39|) with respect to the time, where d'd/dxQpUh and 
xd'^'d / dxQpdxoq are considered to be constant 

-T^^/^ + XTr-^T— -tUt (40) 



A V dxop dxopdxoq J ^^^^ A \dt / ^^^^ 

This assumption is justifled by the fact t hat, according to the classical for- 
mulation of motion of continuum media ( Truesdell ( 1977f )). dd/dx^pUh and 
xd'^'d / dxQpdxQq are smooth functions of t -at least during the period of a 
fluctuation- whereas the fluid deformation varies very rapidly according to 
Eqs. dSTD-dSHl). 



6. Statistical analysis of temperature derivative 

As explained in this section, the Lyapunov analysis of the local deforma- 
tion and some plausible assumptions about the statistics of u and lead 
to determine the distribution function of dd/dr and all its dimensionless 
statistical moments. 

The statistical properties of dd/dr, are here investigated expressin g veloc- 



i ty an d temperature through the following canonical decomposition (jVentsel 



mm) 
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where and 0^ are proper coordinate functions of t and x, and {k = 
1,2, ...) are dimensionless independent stochastic variables which satisfy 

{^k)=0, {^^Q = 6,„ Uj^k) = r^^Jk P, bl»>l, (4') = 0(1) (42) 

where Wijk = 1 for i = j = k, else vuijk=0. It is worth to remark that the 
variables are properly chosen in such a way that they express the mech- 
anism of cascade for both velocity and temperature, through the condition 
\p\ »> 1. 

The dimensionless temperature fluctuation is obtained in terms of ^ 
substituting Eq. (gl]) into Eq. fHO]) 

ij k 

where r = "tAt, d = d 9, whereas Pe = R Pr and R = u\t/v are Peclet 
and Reynolds numbers refereed to the Taylor scale, and Pr = is the 
fluid Prandtl number, therefore Ylij ^ij^i^j l/Pe bk^k arise from con- 
vective term and fluid conduction, respectively. Now, tha nks to the loca l 



isotropy, both u and are two gaussian stochastic variables (jVentsell (119731 ). 



Lehmann (jl999h ). accordingly, satisfy the Lindeberg condition, a very gen 



eral necess a ry an d sufficient condition for satisfying the central limit theorem 



( ^Lehmanru (Il999|)). This condition does not apply to d'dldr = lim A-i^/r. In 

fact, as A-d is the difference between two correlated gaussian variables, its 
PDF could be a non gaussian distribution function. To obtain this PDF, the 
fluctuation dd/dr is first calculated in function of 

f-E^e...i,Ete....S.P.A. (44) 

ij k 

This fiuctuation consists of the contributions L, 5, P and A^, appearing 
into Eq. flH|) : in particular, L is the sum of all linear terms due to the 
fiuid conductivity, S = SijC,iC,j is the sum of all semidefinite bilinear forms 
arising from the convective term, whereas P and N are, respectively, the 
sums of definite positive and negative quadratic forms, which derive from the 
convective term. The quantity L + S tends to a gaussian random variable 
being the sum of statistically orthogonal terms, while P and do not, as 
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they are linear combinations of squares (jMadowl (119401 ). iLehmanru ( 1l999l )). 
Their general expressions are 



P = Po + r]i + 7]i 

(45) 

AT = iVo + Ci - CI 

where Pq and Nq are constants, and r/i, 772, Ci (2 are four different centered 
random gaussian variables which are mutually uncorrelated thanks to the 
hypotheses of fully developed flow and isotropy. Therefore, the longitudinal 
fluctuation of the temperature derivative can be written as 

§ = M+(^2(r/^-l)-^3(C^-l)) (46) 

where ^, 77 and ( are independent centered random variables which exhibit 
gaussian PDF with (^^) = {rj'^) = (C^) = l, and tpi il)2 and ip'i are given quanti- 
ties. Due to the isotropy, the skewness of dd / Or must be equal to zero, thus 
^2 = V's, and 

Furthermore, comparing the terms of Eqs. (jlT]) and (jll]), we obtain that ipi 
and iIj2 are related each other and that their ratio ip = ipi/ip2 depends on the 
Peclet number 

^\ mMj + ^/p^dh/dfi,f) ^ > 

Taking into account the properties (jl2]) of {{il) »> 1, {^t) = and 
in view of Eq. (j48|) , we found 

ilj = — = C^/P~e (49) 

^1 

where C is a proper constant which has to be identified. Hence, the dimen- 
sionless longitudinal temperature derivative is 

^1 + 4^2 ^ ) 
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In order to identify C, observe that Eq. (13U]) is formally similar to the ex 



pressi on of the longitudinal velocity derivative dur/dr obtained in lde Divitiis 



teOlOl ) (see also the appendix) 

dur/dr ^ in + i^u {xiVu^ - 1) - (C' - 1)) 

^{{dUr/dry) VI + 2^2 (1 + ^2) 



(51) 



^u, Vu and (u are independent centered gaussian random variables with (^^) 
(Vl) = iO =1, and 



^„(i?) = ^/— ^V^„(0), 4(0) =0(1), x = x{R) = Oil) (52) 
15v 15 



X 7^ 1 provides a negative skewness of dur/dr, whereas '?/'u(0) ~ 1.075 is 
determined throu gh an approximate estimation of the critical value of R 
(Ide DivitiisI toitj )). Now, when Pr = 1, it is reasonable to assume that the 



ratio between linear and quadratic terms of Eq. fl50|) is equal to that of the 
corresponding terms of Eq. f l5T]) . Accordingly, ip c^l ip^ and this identifies an 
approximate value of C 

C^t^ ^0.141 (53) 

The distribution function of the temperature derivatives is thus expressed 
through the Frobenius-Perron equation 

/ / [p{Opiv)p{C)H€-M^,vX))d^drjdC (54) 

Jri Jc 

where ^r{^, C) is determined with Eq. ( !50l) . 6 is the Dirac delta and p is a 
centered gaussian PDF with standard deviation equal to one. 

Finally, the dimensionless statistical moments of -^r are easily calculated 
considering that ^, rj and ( are independent gaussian variables 

- tBi = 7T-^ ti'l ) *'(?"-'=>(("^ - C^)'> (55) 



(tf2)../2 (1 + 4^.2)»/2 zLj, 



It is worth to remark that, for non- isotropic turbulence or in more complex 
cases with boundary conditions, the stochastic variables C,k could not satisfy 
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the Lindeberg condition, thus ^9 will be not distributed following a Gaussian 
PDF, and Eq. flSOp changes i ts ana lytical form and can incorporate more 
intermittent terms (iLehmann which give the deviation with respect 

to the isotropic turbulence. Hence, the absolute statistical moments of dr 
will be greater than those calculated with Eq. fl55|) . indicating that, in a more 
complex situation than the isotropic turbulence, the intermittency of can 
be significantly stronger. 



7. Self-Similar temperature spectrum 



An ordinary differential equation which describes the spatial evolution of 
fe is now derived from Eg. (|30|1. ado p ting the hypo t hesis of self-similarity 
of Ivon Karman fc LinI (119491 ) . iGeorgd (119881 ) - iGeorgd (Il992l ). and using the 
proposed closure of the Corrsin equation. 

Far from the initial condition, the simultaneous effect of temperature 
and velocity cascade with the fluid conductivity and viscosity, acts keeping 
fe similar in the time. This is the idea of self-preservi ng correlation func- 



tion Av hich was originally introduced by von Karman (see Ivon Karman fc Lin 



(119491 ) and reference t herein) f or what conce r ns the velocity correlation, and 
thereafter adopted by George ( 1988[)- George (|l992h for stud ying the temper- 
ature spectrum. According to Georgel ( 19881 )- George ( 19921 ). the self-similar 
temperature correlation can be scaled with respect to Xrif), thus 



fe = fe{r), where f = r/Xxit) 



Substituting Eq. (1561) into Eq. (1301) . we obtain 



dfe f dXx 
df u dt 



-f dfe 2 
2 dr RPr 



d-r'^ 



^2dfe_ 
f df 



+ 



12 



R Pr 



(56) 



^ ) /.(57) 



Therefore, the boundary problem given by Eqs. ( l30l) and ( 1311) . is here re- 
duced to an ordinary differential equation of the second order in the variable 
r. Equation (|57|) is a non-linear equation whose coefficients vary according 
to Eq. (fTOOD and 



dv^ 
dt 



A ji 



d9^ 

IF 



(5^ 
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Now, if the self-similar ity is assumed, all the c o efScients of Eg. (IFTD do 
not va r y with the time (Ivon Karman fc HowarthI (jl938l )- lvon Karman fc Lin 
fll949h . lGeorg3 fllQSSh - lGeoreel fll992h ). t\ms 



R = const. 



1 dXn 



u dt 



const. 



const 



As At follows Eq. fllOOp . is obtained from the constancy of Xg/Xr 



Xg{t) = Xg{0)^l + 10u t/XliO). 



(59) 



(60) 



Thus, according to Warhaft fc Lumley ( 1978 ) and George f 1988 )- Georgel 
(I1992I ). the microscales Xt, Xg and the rates dO"^ /dt and dv? /dt, depend on 
the initial conditions of temperature and kinetic energy spectra. Taking into 
account Eq. fllOOp . we obtain 



IdXr _^ 
u dt R 



(61) 



Therefore, fg{r) obeys to the following non-linear ordinary differential equa- 
tion 



^ dfg 

— — r + 
R df 



l-fdfg^ 2 
2 df R Pr V df^ 



d^fg ^2dfg_ 

f df 



J^fXr, 
RPr\x,' 



0(62) 



The self-similar solutions are searched over the whole range of f, but for 
the dimensionless distances whose order magnitude exceed R. This corre- 
sponds to assume the self-similarity for all the f requencies of the energy 
spectr u m, with the exception of the lowest ones (jvon Karman fc Howarth 
(1938), von Karman fc Lin ( 19491 )). Accordingly, dfg/dt, can be neglected 
with respect to the other terms 



2 dr \ dr'^ r dr I \ Xg 







The boundary conditions of Eq. (1631) arise from Eqs. (|TT 

/.(O) = 1 



(63) 



(64) 



lim /e(f) = 



(65) 
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For f = 0, Eq. (1^ gives 



a; 



(66) 



To determine fe, Xr/^e must be first specified into Eq. (!63l) . This is de- 
termined using again the self-similarity of / and fe, thus Xr/Xg is calculated 
by substituting Eq. (!60|) into Eq. (!58l) and integrating this latter from t = 
to t 

2 1 



In 



m 

e{o) 



_in(l + 
5 V A J Pr ' 



lOz/ 



whereas the velocity standard deviation is 

uit) 



In 



u{0) 



2 



1 + 



lOz/ 



At(0) 



(67) 



(6^ 



The full self-similarity (mechanical and thermal) occurs when 6 and u are 
proportional each other 

e{t) u{t) 



9{0) m(0) 



(69) 



The value of Xt/Xq satisfying this condition depends on the Prandtl's num- 
ber, and is calculated with Eq. ( 1671) 



Xx 

Accordingly, /g (0) is related to Pr 

d'feiO) 



6 1 
5P^ 



-Pr 



(70) 



(71) 



df'^ 3 

This r esult, in agreement with Corrsin ( JAP 1951 ). George ( 1988 ). Georgel 
(119921), expresses a further link between / and fg only in the case of self- 
similarity. 

Observe that the solutions fe G [0, oo) of Eq. ([63]) with dfg/df{0) = 
and Pr ^ 0, tend to zero as r — )■ oo, thus the boundary condition (l65l) can 
be replaced by the following condition in the origin 

dfg{^) 



df 







(72) 
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Therefore, the boundary problem represented by Eqs. (1^ . (1M|) and (155]) . 
is reduced to the following initial condition problem written in the Cauchy's 
normal form 



— =Fe 
df 




(73) 



the initial condition of which is 

/,(0) = 1, F,(0) = (74) 

In conclusion, the self-similar functions / and fe are calculated as the 
solutions of the ordinary differential system (I102p and (|73l) with the initial 
conditions ffTOal) and (171). 



8. Results and Discussion 

The self-similar temperature and longitudinal velocity correlations are 
here calculated with Eqs. fll02p and (1731) . for several values of R and Pr. 

The case with Pr — )■ (infinitely conductive fluid) is first considered. 
This is a limit case of the differential system fl73|) - fl7^ corresponding to the 
following equation 

^ + -^ = (75) 

df"^ f df 

which does not admit analytical solutions for the boundary conditions fl3ip . 
Conversely, the case Pr — )• oo {u »> x) gives the following result 

lim — — - = — oo (76) 

f^o dr^ 

which expresses the behavior of the temperature correlation near the origin, 
whereas, for f > 0, is obtained solving Eqs. (I751)- (rn|) by quadrature, in 
terms of / 

In _12^f_i^ (77) 

R Jo v^WlO 
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This equation states that, if R is large enough and 



/ \ 2/3 / \ 2/3 

/ - 1 - i-^j , then also ~ 1 - ( £j (78) 
where L0= are length scales proportional to Xt 

Therefore, in case of self-similarity, with Pr — )■ oo, the ratios between the 
scales are 

^ = 1, ^=0 (80) 

When R and Pr change, the ratios between the scales vary depending on 
the combined values of R and Pr, therefore quite different situations occur. 

In order to study the influence of R and Pr on /e, Eqs. f ll02p and ( 173|) are 
numerically solved for different values of R and Pr. The Reynolds number 
is assumed to be R =50, 100 and 300, whereas Pr ranges from 0.001, to 10. 

Figure [3] shows / (dashed hues) and fe (solid line) in such these condi- 
tions. The temperature correlation, related to / by means of the mechanism 
of temperature cascade (see Eq. (l28l) ). is furthermore linked to / by self- 
similarity (1701) . Therefore, for assigned values of R and At, the Corrsin 
microscale decreases with Pr and the curves of fo seem to collapse into a 
single diagram when Pr — )■ oo. On the contrary, small values of Pr, repre- 
senting high thermal conductivity, determine large scales of variations of fg 
and p*. In particular, the case P = 50 is first considered. For Pr =0.001, 
fg exhibits oscillations whose amplitude decreases when f rises. As Pr in- 
creases, the oscillations magnitude diminishes, and for 0.01 < Pr < 0.1 these 
oscillations vanish, being fg > 0, whereas the integral scales and the Corrsin 
length diminish. The case R = 100 differs from the previous one. In fact, 
the higher value of R determines sizable reduction of the oscillations, whereas 
the integral scales of / and fg are greater than the previous ones. Next, for 
R = 300, the integral scales increase again, resulting fg > a monotonic 
function of r for each value of Pr. 

Accordingly, also the triple correlation temperature-velocity varies 
with R and Pr. For R = 50, small values of Pr (0.001) cause large scales 
of variations and sizable oscillations of p*, whereas higher Prandtl numbers 
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produce the loss of these oscillations and a reduction of the length scales and 
of Increasing R{R= 100 and 300), the length scales rise, the oscillations 
disappear, and a reduction of |p*|max is observed. The Prandtl's number 
acts on fg, in such a way that its increment causes a diminishing of the length 
scales and of When Pr = 0.7 a nd 1, the o b tained results agree very 
well with the classical experiments of iMills et all (Il958[ ) which regards the 
turbulence behind heated grid, in the sense that, fg is roughly equal to /, 
whereas p^, and the triple velocity correlation k exhibit the same order of 
magnitude. 

As far as the velocit y correlatio r is, f( r) and k{r), are concerned, these 
agree with the results of Ide Divitiid (1201 ll ). 

The properties of G(k) are the consequence of the temperature correla- 
tions. These spectra, calculated with Eq. fl33l) and flM|) . are depicted in 
Fig. m The variations of Q{k) with R and Pr are quite peculiar. In any 
case, according to Eq. fl28|) . n — )■ 2, as k — > 0. For Pr = 0.001, when R 
ranges from 50 to 300, the temperature spectrum shows essentially two re- 
gions, which correspond to two different scaling laws 0(k) ~ k'^ (see also 
Fig. E]): one near the origin where n ~ 2, and the other one, at higher k, 
where —17/3 < n < —11/3, (value very close to —13/3). Between these re- 
gions, the exponent n varies rapidly at low Reynolds number, whereas when 
R increases, n exhibits more gradual variations. The value of n ^ —13/3 
here c alculated, is in between the exponent proposed by Batchelor et al 
fll959h (-17/3) and the value determined bv iRoeallo elal (Il989h f-ll/iT 
with the numerical simulations. Increasing k, n strongly diminishes, and 
G(k) does not show scaling law. When Pr=0.01, the three curves intersect 
each other for n ~ —5/3, where these have inflection points. The width 
of the interval where —17/3 < n < —11/3 di minishes, in pa r ticula r now 
— 17/3 < n < —13/3, value in agreement with iBatchelor et all ( 119591 ). For 
Pr =0.1, the previous scaling law with n ^ —13/3 vanishes, whereas for R = 
50 and 100, n changes with k, and G(k) does not show a noticeable scaling 
law. When R = 300, the born of a small region in which n has an inflection 
point is observed for n ^ —5/3. For Pr = 0.7 and 1, with R = 300, the width 
of this interval is increased, whereas at Pr = 10, and R = 300, we observe 
two regions: one interval where n has a local minimum with n ~ —5/3, and 
the other one where n exhibits a relative maximum, with n ~ — 1. For larger 
K, n diminishes and the scaling laws disappear. 

Figure m reports also (on the bottom) the spectra T{k) (solid lines) and 
T{k) (dashed lines) which describe the mechanism of kinetic energy and 
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temperature cascade. As these latter do not modify the value of 6 and u, 
Jq T{K,)dK = 0, and r{K,)dK = 0. 

The presence of the scaling law n ~ ~5/3 a grees wit h the theoretical argu- 
ments oflCorrsinl (IJAP 19511 ). lObukhovl (119491 ) (see also lMvdlarski fc Warhaft 
fll998h . lDonzis et all fl201of r and references therein) , where in the limit of high 
R and Pr 



QU) = C,e-^/3 -5/3 



in the inertial-convective range, being 



u 



^1) 



(82) 



and Ce is the so-called Corrsin-Obukhov constant, a quantity of the order of 
the unity. This study identifies Cq by means of the obtained results, analysing 
the quantity 



0(/c)ei/3^5/3^^-i 



(83) 



-here called Corrsin function- in terms of k. Re and Pr. This is calculated 
as Ce = Fc{k,) in the range where Fc{k) is about constant. A different 
Corrsin-Obukhov constant Cei can be also defined with respect to the one 
dimensional spectrum 



dQi 
dKi 



2ki 



where ei(K) = Ceie^^/^egK"^/^ and Cei = 0.3Ce (84) 



Figure E] reports Fc in terms of k for different values of R and Pr. For 
Pr =0.01, at relatively small Reynolds number, the temperature spectrum 
does not follow k~^^^, thus Cg is not defined, wheres at i? = 300 the diagram 
shows a region with a local maximum where the variations of Fq are relatively 
small. This maximum can identify the value of Cg which results to be about 
1.5 [Cgi 0.45). For Pr =0.1, the larger scaling interval, implies a wider 
range in which Fc{k) ~ const, for each Reynolds number, resulting now 
Cg ~ 1.8 {Cgi ~ 0.54). When Pr = 0(1) (0.7 in the figure), the scaling 
law determines that Fc slightly rises with k, ranging from 1 . 4 to 2 
(Cgi ~ 0.42 ^ 0.58), vales comparable with iMvdlarski fc WarhaftI fll998h . 
For Pr = 10, the region with ^ increases, therefore Fq reveals sizable 
variations and Cg can be defined only in a small range of k. 
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The scalinK l aw k ^ is in hne with the theoretical arguments proposed in 
Batchelorl (119591 ). where 



6(k) = Cb\/ - te K, \ 



(85) 



in the viscous-convective range, being Cb = 0(1) the Batchelor's constant. 
The present analysis identifies Cb through the temperature spectra previ- 
ously calculated. To this end, the following quantity 



FbM 



Q{k)k\I — e 



-1 



(86) 



here called Batchelor's function, is considered for several k. Re and Pr. This 
constant is here calculated as Ce = Fb{k) in the region of k where Fb{k) ~ 
const (or at least exhibits a plateau). Also the Batchelor's constant Cbi can 
be defined with respect to the one dimensional spectrum ([8] 



(87) 



0i('ti) = C'si-y/— £0 /^i ^5 where Cbi = 0.5Cb 



In Fig. [71 Fb is represented versus k for different values of Re and Pr. It 
is apparent that Fb ~ const when Pr is high enough. For Pr =10, when R= 
50, 100 and 300, the corresponding values of Cb are about 5, 7 and 8 (that 
is Cbi — 2.5, 3.5 and 4) an d this occurs for 1 < n < 10. These values are in 
quite good agreement with iDonzis et all (120101) farid refer ences therein), and 
are consistent with the experiments of I Grant et all (119681 ) and lOakeyl (119821 ) 
which deal with the temperature spectrum observed in ocean. 



Next, in order to analyse the statistics of the temperature derivative, the 
PDF of dd/dr is calculated with Eqs. (I54p and ( I50p . for different values 
of the parameter = C^Pr R. This PDF is obtained with sequences of 
the variables ^, t] and each generated by a gaussian random numbers 
generator. The distribution function is then calculated through the statistical 
elaboration of these data and Eq. ( l50|) . The corresponding results are in Fig. 
[HI where the PDF is shown in terms of the dimensionless abscissa 



d-djdr 



^{d^/dr^) 



(88) 
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These distribution functions are normalized, in order that their standard 
deviations are equal to the unity. These PDF are even functions of s and 
their tails rise with in such a way that the intermittency of dd/dr increases 
with if), according to Eq. (|55|1 . The PDFs shown in the figure are calculated 
for ip = 0.25, 0.5, 1., 10., oo, in particular for ip = 1Q and oo, the curves are 
about overlapped. In order to study the influence of if) on this intermittency 
and on the statistics of s, the flatness H^^ and the hyperflatness H^, defined 
as 

are also calculated. These are depicted in Fig. |8]in function of ip- For ip =0, 
the PDF is gaussian which corresponds to H4 = 3 and = 15. Increasing 
ip, the non-linear fluctuations due to r] and (, determine an increment of 
and Hq that tend to the limits H4 = 9 and Hq = 225 whe n ip — )■ 00. 



These results are compared with the experiments of ISreenivasan et al 



( 119801 ) which in turn give the flatness of dd/dr. The value of C identified 
with this comparison of is C ~ 0.135 against the value C ~ 0.141 here 
calculated in the proper section. 



9. Conclusions 



The finite scale Lyapunov theory is adopted to study the temperature 
fluctuations in homogeneous isotropic turbulence. This analysis leads to the 
closure of the Corrsin equation and provides the statistics of the temperature 
fluctuations. 

The r esults, whi c h rep resent a further applic ation of the analysis pre- 
sented in Ide DivitiisI ( 120101 ) and Ide DivitiisI ( 1201 ll ). are here obtained in the 
case of self-similar velocity and temperature fluctuations, and can be so sum- 
marized: 



1. The energy equation, written using the referential coordinates and the 
Lyapunov analysis of the local deformation, allows to factorize the tem- 
perature fluctuation and to express it in Lyapunov exponential form of 
the local deformation. 

2. The finite scale Lyapunov analysis provides an explanation of the phys- 
ical mechanism of temperature cascade and gives the closure of the 
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Corrsin equation. This is a non- diffusive closure expressing G in terms 
of / and dfe/dr. 

3. This closure provides a mechanism of cascade which generates tem- 
perature spectra with different scaling laws, depending on R and Pr. 
In particular, for proper values of R and Pr, these spectra satisfy the 
Corrsin-Obukhov and Batchelor scaling laws in opportune regions of 
wave-numbers. 

4. The Corrsin-Obukhov and Batchelor constants here identified with the 
proposed theory, agree with the different source from the literature. 

5. The PDF of dr and the corresponding dimensionless moments are deter- 
mined through a canonical decomposition of velocity and temperature 
in terms of random variables which describe the mechanism of cascade. 
This is a non-Gaussian PDF whose intermittency increases with R and 
Pr, in agreement with the experiments of the literature. 
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11. Appendix 

For sake of convenience, this sec tion reports themain results of the finit e 
scale Lyapunov analysis obtained by lde Divitiid ( l2010l ) and lde DivitiisI ( l201ll ). 
which deal with the homogeneous isotropic turbulence. 



11.1. Closure of the von Karman-Howarth equation 

For fully developed isotropic homogeneous turbulence, the pair correlation 
function f of the longitudinal veloc i ty Ur, satisfies the von Karman-Howarth 



equation (jvon Karman &: HowarthI (119381 )) 



dt 



K{r) 



+ 2u 



+ 



lOz/ 

A ji 



/ 



(90) 
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the boundary conditions of which are 



/(O) = 1, 
hm /(r) = 



(91) 



where r is the separation distance, At = a/— l//"(0) is the Taylor scale, and 
u is the standar d deviation of Ur, which s a tisfies the equation of the turbulent 
kinetic energy ( von Karman fc Howarth ( 19381 )) 



lOz/ 

A ji 



■u 



(92) 



This equation, obtained putting r = into Eg. ([90|). gives t he ra te of kinetic 
energy in function of u and At (Ivon Karman &: HowarthI (119381 ) . iBatchelor 
( Il953l )). The function -ft'(r), related to the triple velocity correlation func- 
tion, represents the effect of the inertia forces and expresses the mechanism 
of energy cascade. Thus, the von Karman- Howarth equation provides the re- 
lationship between the statistical moments ((Am^)^) and ((Am^)^) in function 
of r, where Aur is the longitudinal ve locity differ e nce. 

The Lyapunov theory proposed in de Divitiii ( 2010 ) leads to the closure 
of the von Karman-Howarth equation, and expresses K{r) in terms of / and 
df /dr 



K{r) 



u 



l-fdf 
2 dr 



(93) 



where K{0) = represents the conditio n that K does not modi fy the fluid 
kinetic energy (jvon Karman &: HowarthI (119381 ) , iBatchelorl (119531 ) ) . 



11.2. Statistics of the longitudinal velocity difference 

Here, the results of Ide Divitiis (l2010l ). dealing with the statistics of the 
longitudinal component of velocity difference Aur are recalled. 

There, Aur is represented in terms of centered random variables 



AUr ^ U + {xiVu^ - 1) - (Ca' - 1)) 



(94) 
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where ipu is a function of r and of the Taylor scale Reynolds number 



R 



15V15 



(95) 



"ipuo = il> v.(R,0), with tjjviO ) = 1.075, and x 7^ 1 provides a nonzero skewness 
of Aur f ide DivitiisI fl2010h l 

Equation flM|) . arising from statistical considerations about the Fourier- 
transformed Navier-Stokes equations, expresses the internal structure of the 
fully developed isotropic turbulence, where and Cu are independent 

centered random variables which exhibit the gaussian distribution functions 
p{Cu), vijlu) and p(Cu) whose standard deviation is equal to the unity. 



11.3. Self-Similarity in homogeneous isotropic turbulence 



Now, the results of Ide DivitiisI (120111 ) are briefly summarized, for what 
concerns the self-similarity of homogeneous isotropic turbulence. These re- 
sults are based on the i dea of self-preserving corr e lation function which was 
origin ally proposed by Ivon Karman fc Howarthl ( ll938l )- lvon Karman fc Lin 
( 119491 ): far from the initial condition, the combined effects of energy cas- 
cade and viscosity act keeping the velocity correlation function / and the 
energy spectrum E{k), similar in the time for large values of wavelengths, 
in particular in the inertial sub-range. The condition of self-preserving /, 
applied to Eq. (!90|) . leads to the following ordinary differential equation 



2 dr R 



f dr J R 







(96) 



whose boundary conditions are from Eqs. 
(Il938h ) 



(19T|) ( Ivon Karman fc Howarth 



/(O) 



(97) 



lim /(f) = 

f—KX) 

Into Eq. (I96p . / = f{r), where f is the dimensionless variable f 
which depends upon r and t, therefore 



(if 2 



(0) = -1 



(98) 
r/Xrit) 

(99) 
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Observe that Eq. fl96l) gives the self-similar / and E in the inertial range 
and for k — )■ oo (small r). Although Eq. does not describes the energy 
spectrum near the origin, E{k,) satisfies the continuity equation, being E ^ k,"^ 
for kXt « 1. 

This similarity and the equation of the turbulent kinetic energy lead to 
the expressions of u and At 



Xrit) = At(0)7i + 10i/ t/X^O), u{t) = —= 

V 1 



«(0) 



+ lOi/ t/A|(0) 



(100) 



As the solutions f & [0, oo) with df /df{0) = tend to zero when r 
oo, the boundary condition ( 198|) can be replaced by the following condition 
in the origin 

*P = (101) 

dr 

Therefore, the boundary problem represented by Eqs. ( I96l) . (I97jl and 
(l98l) . corresponds to the following initial condition problem written in the 
Cauchy's normal form 

- F 

df 

(102) 



dF ^ 1 1-f 4\ 

the initial condition of which is 

/(0) = 1, F(0) = (103) 
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Figure 3: Correlation functions for Pr= 10"'^, 10~^, 0.1, 1.0 and 10, at different Reynolds 
numbers. Top: velocity correlation / (dashed line) and temperature correlation fg (solid 
lines). Bottom: triple velocity correlation k (dashed line) and triple velocity-temperature 
correlation (solid lines) 
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Figure 4: Spectra for Pr= 10^'', 10^^, 0.1, 1.0 and 10, at different Reynolds numbers. 
Top: kinetic energy spectrum E{k) (dashed line) and temperature spectrum G(k) (solid 
lines). Bottom: velocity transfer function T(k) (dashed line) and temperature transfer 
function r(«;) (solid line) 
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Figure 5: Scaling exponent of the temperature spectrum calculated for Re = 50, 100 and 
300, at different values of the Prandtl's number. 
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Figure 6: Corrsin function for R= 50, 100 and 300, at different values of Prandtl's numbers. 
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Figure 9: Dimensionless statistical moments, H4 and Hq of d'&fdr in function of the 
parameter tjj. 
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